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Statistical mechanics of one-dimensional s = | anisotropic XY 
model in transverse field with Dzyaloshinskii-Moriya interaction 

For ID s = I anisotropic XY model in transverse field with 
Dzyaloshinskii-Moriya interaction using Jordan- Wigner trans- 
formation the thermodynamical functions, static spin correla- 
tion functions, transverse dynamical spin correlation function 
and connected with it transverse dynamical susceptibility have 
been obtained. It has been shown that Dzyaloshinskii-Moriya 
interaction essentially influences the calculated quantities. 
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1 Introduction 



In 1961 E.Lieb, T.Schultz and D.Mattis in ref. [|T| pointed out one type of 
exactly solvable models of statistical mechanics that is so called ID s = | 
XY models. Rewriting the Hamiltonian of such chain 

^ = E [(1 + 7)^ •+! + (1 - 7)44+i] ' -1 < 7 < 1, (1) 

j 



zSjmsZi, Ci, P,l = x,y, z + cyclic permutations (2) 



with the help of the raising and lowering operators Sj = s| ± iSj in the form 



(3) 



they noted that the difficulty of diagonalization of the obtained quadratic 
in operators s~^, s~ form (3) is connected with the commutation rules that 



these operators satisfy, namely, 



(1 - 2s+s^). Really, they 



are similar to Fermi-type commutation rules for operators at the same site 
and to Bose-type commutation rules for operators attached to different sites 



0; 



0, j 7^ m. 



(4) 



That is why one should perform at first Jordan- Wigner transformation (see, 
besides ref. [0], also refs.[2-4]) 



Cl 



1 ' '^j 

+ „+ 
1 , t-j 



sjPj-i = Pj-is-, j = 2,...,iV, 



AT, 



(5) 



where Jordan- Wigner factor is denoted by Pj = nn=i("~2s^). The intro- 
duced operators really obey Fermi commutation rules. From (5) it follows 
that 



„+p2 - 



(6) 



since Pj = nn=i(~2s^)^ = nn=i4(s^)^ = 1, and the commutation rules at 
the same site remain of Fermi-type. Consider then a product of c-operators 
at different sites 



n— 1 m—1 m—l 

4 n(-24) n(-24)^; 

p=l j=l j=n 



4 n (-2^)^;;^' 



(7) 
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putting here for definiteness n < m. Since sf{—2Sj) = ±sf and {—2Sj)sf — 
=F sf, and consequently 

m— 1 m—l 

Cm4 -s^Ii (-2^1)4 = -4 n (-2s|)s^, (8) 

j=n j=n 

one gets c+c^ = -CmCn- Similarly one finds that c+c+ = -c;t;c+, c„c^ = 
= — CmCn- Thus the introduced in (5) operators are Fermi-type operators 

[cj, 4} = 6ji, {c+, C;+} = {cj, ci} = 0. (9) 

Since Pj = 1, Pj = exp(±«7rX;n=i ^n^n) (because exp ±^7^X]n=l(| + ^n) = 
= nn=i(~2s^)), Sj'sJ — cfcj, it is easy to write the inverse to (5) transfor- 
mation 

i-i j-i 
Si=ci, = Cjexp(±z7r^c+c„) = exp(±m^c+c„)cj, j = 2, . . . , A^, 

n=l n=l 

i-1 i-i 
s| = c]^, = c+exp(±z7r^c+c„) = exp(±z7r^c+c„)c+, j = 2,...,N. 

n=l n=l 

(10) 

Returning to the Hamiltonian (3) one notes that the products of two 
Pauli operators at neighbouring sites transform into such products of Fermi 
operators 

'^/(~24)'4+i = ^t^t+i^ 
s+{-2s^)sJ^, = sp-^,, 

Usually bearing in mind the study of thermodynamical properties of the 
system that requires the performance of thermodynamical limit N oo, 
the periodic boundary conditions are implied 

s%+i = Si, a^x,y,z. (12) 

In connection with this in (3) there are products of the form 

N-l 

— ^N^l — ^l^N — ^l^^N n (~2Sp) = cfcj^P, P = Pn, 

p=l 

~ Clcj^P, "S^v^iV+l ~ ~C^C]\[P, Sj\fSj^_^_i = —CiCjsiP. (13) 





-1 






-1 








s 


CjCj^l 




s 
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Gathering the obtained terms one finds that for the ring 
H^H- + BP+ = H+P+ + R-P-. 

Here 



N 



(14) 



(15) 



the difference between and H is only in the imphed boundary condi- 
tions: for they are antiperiodic 



^3 



and for H they are periodic 



(16) 

(17) 



B = 



H- 



{^Cj^Ci + c^fCij + h.c. is the boundary term; P^ 



(l±P)/2 are orthogonal projectors (P++P- = 1, {P^f = P^, P^P^ = 0), 
besides this [H^,P] = [H^,P^] = 0. For the open chain with free ends (then 
in the sum in (1) the range of the summation index is j = 1, . . . , A'^ — 1) the 
Hamiltonian after fermionization has the form 



N-l 



H^IY. [(744+1 + c+c.+i) + h.c. 



(18) 



Formulae (14), (15) or (18) realize the reformulation of the initial Hamil- 
tonian (1) in terms of fermions. They are the starting point for further study 
of statistical mechanics of models hke (1). Besides it appears [5,6] that for 
calculation of free energy 



/ 



-— lim 

/3 iV^oo 



— Spexp{-f3H) 



(19) 



or static spin correlation functions 



s]: >= Ihn {Sp [eM-PH)s'^' ■ ■ ■ ^1"] /^pexp(-/3i/)} (20) 



the boundary term may be omitted and hence one has to consider a system 
of free fermions. It is more difficult to calculate the dynamical correlation 
functions. Really, 

Sj{t) = exp{iHt)sj exp{-iHt) = 

= exp{iH+t)s^j exp(-2i7+t)P+ + exp{iH-t)s^j exp{-iH-t)P~ = 

= P+ exp{iH+t)s'j exY>{-iH+t) + P" exp{iH-t)s^j exp(-«i/"t) (21) 



5 



(owing to the following relation that is valid for arbitrary function of H = 

H+P+ + H-P-. fiH) ^ EZo ^{H^P^ + H-p-r = En=o ^ X 
x[{H+YP+ + {H-YP-]=f{H+)P+ + f{H-)p- = P+f{H+) + p-f{H-)) 
in contrast to 

sY^i) = 

= exp{tH^t) sf" exp{-tH^t)P' + exp(«i/"t) sf^ exp{-tHH)P+ = 
= P+exp{iH+t) sf' exp{-iH-t) + p-exp{iH-t) sf^ exp{-iH+t). (22) 

In accordance with (21) the pair transverse correlation function in the ther- 
modynamical limit can be written as 

^ _ Sp [exp{-pH-) exp{iH-t)s'j cxp{-iH-t)s'j^^^ 
< sAt)s,,^ >— [26) 



^■+" Sp exp {-/3H- 

and hence may be calculated with c-cyclic Hamiltonian. Whereas the pair 
longitudinal correlation function in accordance with (22) in the thermody- 
namical limit can be written as 

_ Sp [exp{-pH-)exp{iH-t)s^exp{-tH-t)0-{t)s^^^ 
<sj{t)sj^^>- Spexp{-PH-) 

(24) 

where O (t) = exp{iH t) exp {—i{H + B)t). The calculation with c-cyclic 
Hamiltonian that neglects the boundary term B yields the approximate 

result that, in particular, is incorrect in the limit of Ising model (7 = 1) 
(see, for instance, [7]). It is interesting to note that the calculation of the 
four-spin correlation function in the thermodynamical limit involves only 
c-cyclic Hamiltonian 

Sp exp{—PH~) exp{iH~t)s^^Sj^ exp{—iH~t)Sj^s 



Sp exp {-PH-) 



(25) 



Thus here as in the case (23) one comes to calculation of the dynamical cor- 
relation functions of the system of non- interacting fermions (see [8]). The 
calculation of the pair longitudinal correlation function, in spite of a great 
number of papers dealing with this problem, remains an open point of sta- 
tistical mechanics of ID s = | XY models. Among other interesting and 
principal questions of the theory of ID s = | XY models one may mention 
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the investigation of nonequilibrium properties of such models (see, for ex- 
ample, [9]) and the examination of the properties of disordered versions of 
such models (see, for example, [10]). 

It is necessary to stress the essential features of the present consideration: 

• the dimension of space D=l; 

• the value of spin s = |; 

• interactions occur only between neighbouring spins (otherwise the Ha- 
miltonian will contain the terms that are the products of more than 
two Fermi operators); 

• only X and y components of spins interact and the field that may be 
included should be transverse (the interaction of z components, for 
instance, leads to the appearance in the Hamiltonian of the terms that 
are the products of four Fermi operators). 

In connection with this it is easy to point out the model that has more 
general than in (1) form of interspin interaction, and that still allows the 
described consideration. Really, considering the additional terms in the 
Hamiltonian that have form J2j ( J^^sJsJ_,_i + J^^s^sJ_,_^ j one notes that after 
fermionization they do not change the form of the Hamiltonian (14), (15) or 
(18), and lead only to changes in the values of constants. The Hamiltonian 
of the generalized ID s = | anizotropic XY model in transverse field that 
as a matter of fact will be studied in the present paper is given by 



Before starting the examination of this model it is worthwhile to mention 
its possible physical application [11]. For this purpose let's perform the 
transformation of rotation around axis z over an angle a 



Then rewritting at first new terms in sum in the Hamiltonian (26) in the 
form 



3 



3 



Sj — s| COS a + sin a, 
Sj — s| COS a — s| sin a. 



s| = — s| sin a + cos a, 
= s| sin a + s| cos a, 





(28) 
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taking into account that the terms (^SjS^_^_i — s|sj_|_i j are invariant under 
transformation (27) and 

jxy _|_ jyx 

J^^ cos^ a + Z sin 2a + J^^ sin^ a ] sjsj+i + 



+ 



jyy _ j^x j'^y + jy^ \ / r V 
sm 2a H cos 2a I isj Sj_^i + SjS 



7.y 7,x 
i+i 



+ 



+ J^^ sin^ a - 



jxy^jyx _ 



sin 2a + J^^ cos^ aj s^js'j_^i , 



(29) 



and choosing the parameter of transformation a from the condition ( J^^+ 
+jy^) cos2a—{J^^ — Jyy) sin 20; = 0, one will have 



(30) 



where 



jxy _|_ jyx 

jx ^ jxx ^Qg2 ^ _^ ^ g-j^ + J^'?' sin^ a, 

jxy I jyx 

jy = sin^ a sin 2a + J^^ cos^ a, 

T2:y _ Tyx pj/ _|_ jyx 

D=- tB.n2a^ ^ . (31) 

2 Jxx _ Jot ^ ^ 



In the term that is proportional to D one easily recognizes z component of 
the vector [sj x Sj+i] that is the so called Dzyaloshinskii-Moriya interaction. 
It was at first introduced phenomenologically by I.E. Dzyaloshinskii [12] and 
then derived by T.Moriya [13] by extending Anderson's theory of superex- 
change interactions [14] to include spin-orbital coupling (see, for example, 
rcf.[15]). The model with relativistic Dzyaloshinskii-Moriya interaction to- 
gether with ANNNI model are widely used in microscopic theory of crystals 
with incommensurate phase [16,17]. In the classical case Dzyaloshinskii- 
Moriya interaction may lead to the appearance of the spiral spin structure. 
The possibility of the appearance of spiral structure in quantum case has 
been studied in ref.[ll] where for this purpose pair static spin correlation 
functions have been estimated. 

Except the mentioned paper [11] the problem of statistical mechanics of 
ID s — \ XY type model with the Hamiltonian (26) or (30) as far as the 
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authors know was not considered yet ^. At the present paper an attempt to 
fill up this gap by the generalization for this case of the well-known scheme 
of consideration of ID s = | XY model has been made. In section 2 the 
transformation of the Hamiltonian to the initial form for further examina- 
tion of statistical properties is presented. In section 3 the thermodynamical 
properties of the model are considered, and in section 4 it is shown how to 
calculate the static spin correlation functions in this model. The dynamics 
of transverse spin correlations and the transverse dynamical susceptibility 
are studied in section 5. The conclusions form section 6. 



2 Transformation of the Hamiltonian 

In the spirit of above described approach the Hamiltionian of the model 
(26) at first should be rewritten with the help of the raising and lowering 
operators in the form that is similar to (3) 

^f = *^g(»K-2) + 

N 

J=l 

J++ ^ ^jxx _ jyy _ ^^jxy ^ jyx^^ ^4 _ (J"")*, 

j+- ^ ^jxx ^ jyy ^ ^^jxy _ jyx^^ ^4 _ ( J"+)*; (32) 

here the periodic boundary conditions (12) are imposed. The Hamiltonian 
of the model (32) after Jordan- Wigner transformation (5), (10) will have the 
form that is similar to (14), (15) 



H = H+P+ + H-P-, 

3=1 ^ j=l 



—J ~^CjCj^i — J CjCj+i) ; (33) 



besides H~ is c-cyclic and is c-anticyclic quadratic forms in Fermi oper- 
ators. After Fourier transformation 



^In ref.[18] on the base of the model with Hamihonian (26) the problem about the va- 
lidity of the Bose commutation rules approximation for spin operators has been examined. 
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E 



e^^^c^ c = 



(34) 



with = 2Tm/N for i/" and = = 27r(n + l/2)/iV for 

n = -7V/2, -Ar/2 + 1, . . . , A^/2 - 1 (for even N), n = -(iV - l)/2, -(AT - 

l)/2 + 1, . . . , (AT - l)/2 (for odd N) can be rewritten in the form 



-- + e«c+c« - z sin «; ( J++c+c+^ + J c«,c_«,) 



cos AC, ^ = 



jsj/ _ jyx 



sin K (35) 



2 ' 2 

(here the following relations 

Y e"'" c+cl^ = E c+cl^, 5Z e'" c«c_„ = z sin k c«c_« (36) 

were used). The diagonalization of the quadratic forms is finished up by the 
Bogolyubov transformation 

A, = y,y*_^-x,x*_^^0. (37) 
/3-operators remain of Fermi type if 

Kr+ii/.r = i, ^+^=0. (38) 

y-K 

The transformed Hamiltonian contains the operator terms proportional only 
to if 

(39) 



The condition (39) and the second condition in (38) yield 



2IJ++I siuK 



6^+^)2 + 4|J++|2sin^/«. 



(40) 



Taking into account the first condition in (38) one finds that for lower sign 
in (40) 



X 



« = 2z| J++I suiKeK.^{-iargJ++)/\l2£^{£^ - e^+^), 

yn = ^{S,-e^^^)/28,; 



(41) 
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besides = 1, = ^ + For upper sign in (40) 



= -2«| J++I sinre Q-x.Y>{-iargJ^^)/\J 2£^{8i^ + ek^^), 

= V(£« + el+V2^K; (42) 

besides = 1, E,^ — e^~^ — 8,^. Thus in a result of Bogolyubov transforma- 
tion (37) one gets 

H^ = Y.E. - 1/2) , (3^) = S,^, {P,, /3J = = 0. 

K 

(43) 

It is important to note that in contrast to anisotropic XY model because 
of inequality J^^ ^ J^^ one has ^ E_^. This is connected with the 
absence of symmetry with respect to spatial inversion. Really, the Hamil- 
tonian of the model (26) H{Q, J^^, J^^, J^^, J^^) under the action of spatial 
inversion, that leads to change of indexes j to —j oi N —j, j + 1 to N —j — 1, 
transforms into H{Q, J^^, jy"^, J'^y, jyy). 

In fig. 1 it is shown how the presence of Dzyaloshinskii-Moriya interaction 
influences the dependence of E^^ = ej,^^ + S,^ on k in dc Gcnncs model (1 : 
D = 0, n = 0, V : D = J^^, = 0; 2:D = 0, = J^^, 
2' : D = J^^, Q = J^^). In fig. 2 the same is depicted for the case of 
isotropic XY model. 

It is worthwhile to note that the spectrum of elementary excitations in 
the model under consideration as it follows from the expression for ground 
state energy (46) is given by \E^\. 



3 Thermodynamics 

For investigation of thermodynamical properties of the model in question 
let's calculate the free energy per site in the hmit — > oo. In accordance 
with refs.[5,6] for such calculation one can use c-cyclic Hamiltonian and thus 



f = lim 



N 



hiSp exp {—(3H 



(44) 



The diagonalized quadratic in Fermi operators form H involved in (44) has 
the form (43), and owing to this one easily obtains the desired result 



N 



In Sp exp 
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= -- lim 
= -\ lim 

1 1 



^In I TT2cosh^^ I =—4 lim In (2 cosh 



/3 27r i-TT 



In I 2 cosh 



(45) 



Knowing the free energy (45) one finds the energy of the ground state 



= lim / = -— / dK\E^\, 



(46) 



the entropy 
dl3 

the specific heat 



1 r-^ 

~ 2TrJ-n 



c?«; In 2 cosh 



1 

2t: J-TT 



dK^tanh^, (47) 



^813 27r 



UK — — 



cosh 



the transverse magnetization 



N 



< 



N 



on 



1 r , dE I3E, 
-— / an-—- tanh— — . 



(48) 



(49) 



the static transverse susceptibility 



Xzz — 



1 

47r J-n 



do? 



tanh 



Stt 



/TT 
-TT 



cosh 



-2 



.(50) 



In order to illustrate the infiuence of the additional interactions on ther- 
modynamical properties let's present the results of numerical calculations of 
the specific heat (48) as a function of temperature (figs. 3,4 (1: D = 0, f2 = 
0, V : D = J^^, 1] = 0; 2: D = 0, = J^^, 2' : D = J^^, n = J^^)) 
and of transverse magnetization (49) as a function of transverse field (figs. 
5,6 (1: D = 0, P ^ 1' : D - J^^, /3 = 2: 

D = 0, /3 = 1000/ J^^, 2' : D = J^^, p = 1000/J^^)) and of tem- 
perature (figs. 7,8 (1: D = 0, n = 0, 1' : D = J^^, Q = 0; 2: 
D = 0, n = J^^, 2' : D = J^^, n = J^^)) for de Gennes model and for 
isotropic XY model in the presence of Dzyaloshinskii-Moriya interaction. 
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4 Static spin correlation functions 



For the investigation of the spin stucture in the model under consideration 

let's introduce the static spin correlation functions. Due to the possibility 
of their calculation with the help of c-cyclic Hamiltonian the initial formula 
for their evaluation can be rewritten in the form 

< 4^ . . . sj; >= ^im {Sp [exp {-(3H-)sf^ . . . s^^] / Sp exp {-(3H-)} . 

(51) 

Let's introduce then (/9-operators that owing to (37), (41) are the following 
linear combinations of /3-operators 

^ 4 ± 9 = ± , 

K 

\% = -^e^^K^^ ± Vn), = ^e-^[x«exp (2zarg J++) ± y^], (52) 
and in terms of which the spin operators can be presented as 

93-operators obey the following commutation relations 

{(^+ (^7} = 0, {ift, (^+} = 25,,-, {(^-, = -25,,-, (54) 

besides {ipj'ipjy = 1 and 

^t^ ^^^7] = '^^ijVf^ [ftvT, ^t^j] = 0- (55) 

That is why the calculation of static spin correlation functions after substi- 
tution of (53) into (51) and exploiting of (54), (55) reduces to application of 
Wick-Bloch-de Dominicis theorem. The theorem states that the mean value 
of the product of even number of (p operators with the Hamiltonian H~ (43) 
is equal to the sum of all possible full systems of contractions of this product; 
if the number of (/? operators in the product is odd the mean value of the 
product is equal to zero. The full system of contractions of the product of 
even number of Fermi-type operators forms so called Pfaffian the square of 
which is equal to the determinant of antisymmetric matrix costructed in a 
certain way from elementary contractions [19,20]. 
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Thus let's consider the calculation of elementary contractions. One has 

Kl,K2 

(56) 

(here evident relations < P^-^f3i^2 >~< /3_«^j/9_«^2 >= were used). Since 



where = 1/(1 + e^^"^) and in accordance with (52) 
1 



At 



,+ 



2 1 J"*""*" I sin ( argi J"*""*" ) sin k 



one has 



Similarly one finds that 



N 



N 



and 



(57) 



(58) 



< >= ^ E e^'^" (1 + ^«) - I E + A- (59) 



< ^Jcpj^^ >= E e'"" (1 - -5.) + I E«i^ - -5,)/, (60) 



1 



.{+) 



N 



cos (fi;n)-| sin [Kn)Ci^ 



A, (61) 



1 



Ee 



' + J + — E 



N 



cos (ren) h sin (ren) C^^ 



_ 2| J++ 1 cos (arg' J++) sin K 
(-^K = ■ 



(62) 
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The essential simplifications in expressions for contractions (59)- (62) take 
place in the case of model (30), that is when J^^ = —J^^ = D. Then 

< ^t'fj+n ^n,o + H sin (ku) tanh ^ = E{n), 

< ^J^J+n >= -E{n), 
1 PE 



< 'Pi'pt+n >= -G{-n), 



(63) 



where cos = e|+^/£^K, smipK= ^J^^smn/E^. 

Let's return to the evaluation of equal-time spin correlation functions 
and consider, for instance, < >• For this correlation function with 

the utilization of (53)- (55) one derives 



1 



< + ■ ■ ■ ^j+n-lVj+n-lVj+n 



(64) 



and after exploiting Wick-Bloch-de Dominicis theorem in r.h.s. of (64) for 
its square one gets the following expression 



4 < sy^^^ > 



-E{1) .. 


.-E{n 


-1) 


-G(-l) 


-G(-2) .. 


■ -G{- 


n) 


E{1) 


.-E\n 


-2) 


-G(0) 


-G(-l) .. 


.-G{-n 


+ 1 


m E{1) ■■ 


.-eIu 


-3) 


-G(l) 


-G(0) .. 


.-G{-n 


+ 2 


E{n-1) E{n-2) . 







-G{n-2) 


-G{n-3).. 


■ -G{- 


1) 


G(-l) G(0) .. 


. G{n- 


-2) 





E{1) .. 


. E{n- 


1) 


G{-2) G{-1) . 


. G{n- 


-3) 


-^(1) 





. E{n- 


2) 


G(-3) G{-2) . 


. G{n- 


-4) 


-Ei2) 


-E{1) .. 


. E{n- 


3) 


G{-n) G{-n + l). 


■ G{- 


1) 


-E{n - 1) 


-E{n-2).. 








In a similar way for other pair spin correlators one obtains 



(65) 
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-^(1) . 

E{1) 

E{2) E{1) . 

E{n-1) E{n-2) . 

E{n) E{n - 1) . 

G(-l) G{0) 

G{-2) G{-1) 



-E{n-1) -E{-n) -G{-1) . 
..-E{n-2)-E{n-l) -G{0) . 
..-E{n-3)-E{n-2) -G{1) . 

-^(1) -G{n-2). 

.. E{1) -Gin-l). 

.. G{n-2) G{n-1) 

.. Gin -3) Gin -2) -E{1) . 



-G{- 
-Gi- 
-G{- 



-n+1) 
-n + 2) 
-n + 3) 



G{-n+l)G{-n + 2)... G{0) 



G{1) -E{n-2). 



-G(0) 

E{n-2) 
E\n - 3) 





(66) 









:0; 


( 








2 







-E{1) 


-G{1) -G{0) 


-G(-l) . 


..-G(-n+ 1) 


E{1) 





-G{2) -G(l) 


-G(0) . 


..-G{-n + 2) 


m 


E{1) 


-G{3) -G{2) 


. 


..-Gl-n + 3) 


G{1) 


G{2) . : . 


E{1) 


E{2) . 


. . E{n) 


G{0) 


G(l) 


-E{1) 


E{1) . 


.. E{n-l) 


G{-1) 


G(0) 


-E{2) -E{1) 





.. E{n- 2) 


G{-2) 


G{-1) 


-E{3) -E{2) 


-^(1) . 


.. E{n- 3) 



G{-n + 1) G{-n + 2)... -E{n) -E{n - 1) -E{n - 2) 



(68) 



4 < > 



l2 
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-^(1) 


E{1) 





m 


E{1) 


E{n - 1) 


E{n - 2) 


G(l) 


G{2) 


G(0) 


m 


G(-l) 


G{0) 



G{-n + 2)G{-n + 3).. 



-E{n - 


-1) 


-G(l) 


—E[n - 
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(69) 
(70) 
(71) 
(72) 

(73) 



-G(n- 1)... 
E{1) ... 


-E(l) ... 
-E{n - 2). . . 



4 < > 



< s^s^ >- 



0, 



^(n) G'(O) G{n) 
-E{n) G(-n) ^(0) 
-^(0) -G(-n) -£;(n) 



-G'(n) -G(0) £;(n) 

In figs. 9-12 the temperature dependences (1: -D = 0, = 0, \' : D = 
J^^, ^] = 0; 2: D = 0, = J"="^, 2' : D = J^^, 1] = J^^) and in figs. 
13-16 the dependences on transverse field (1: D = 0, (3 = 1/J^^, 1' : D = 
J^^, /3 = 1/J^^ 2 : = 0, 13 = 10/J^^, 2' : D = J^^, /5 = 10/J^^) 
for some pair spin static correlation functions are depicted. It is necessary to 
underline the peculiarities caused by the presence of Dzyaloshinskii-Moriya 
interaction. First, only < >, < >,< >)< -^j-^j+n > 

equal to zero, but not < > and < sjs|^„ >. The last two correlators 

tend to zero when J^^ = J^^ = 0. In this case E{n) = for n ^ and 
hence (66) and (68) may be rewritten as determinants of matrices with 
only non-zero rectangle (but not square) submatrices on their diagonals; 
such determinants are equal to zero. Second, the dependence of pair static 
correlation functions on n is nonmonotonic (in accordance with ref.[ll] this 
fact indicates the appearance in the system of the incommensurate spiral 
spin structure). 
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5 Dynamics of transverse spin correlations 
and dynamical transverse susceptibility 

Let's consider the dynamics of transverse spin correlations calculating for 
this purpose the transverse time-dependent (dynamical) pair spin correlation 
function < Sj{t)sj^^ >. Due to the possibility of exploiting for its calculation 
c-cyclic Hamiltonian H~ (43) the evaluation of this correlation function in 
accordance with (53) and (52) reduces to estimation of dynamical correlation 
functions density-density for the system of non-interacting fermions 

4 < s|(t)s,V >= 

= E < [K.^f^l it) + /^>c./5-i it)] [K^X. it) - )"7«./3-«. it) 



X 



X 



['^j+nA3'^«3 + ^~j+n,K3(^-l^3 ^j+n,K4.(^Ki ^^j+n,K4(^-K4 



>= 



EL \+ \" /,+ /;- ^ R+ R+ R R ^ AEi^l+Ei^2)* ^ 

[ ^jKl^jK2^*'j+n,K3^*'j+n,K4^ h' Klh' K2l^ - 1^31^ - l^i ^ 



K-4,K2,K3,K4 



_\+ ,,- ,,+ \- ^ R+ R R R+ p*(-^Ki--^-«2)* — 

^jKl^^jK2t^j+n,K3'^j+n,K4 ^ A-'kiA^-K2/^-K3/^K4 ^ 

_,,+ A" A+ < B B > p*(^«2--^-''i)* + 

-^/y+ A" /y+ At < /? /?+ /? /?+ > p<^''2-'^-'=i)* _ 

-/^;«4/^7.2^;+n,.3'^7+n,M < /i^-M/i^-«2/i^«+3/i^.+ > e-<^-'^i+^--)*] . (74) 

In r.h.s. of (74) only non-zero averages of /3-operators are written down and 
the following relations 



Ptit) = /3+ exp {lE^t), f3^it) = /3« exp (-lE^t) 



(75) 



were used. The averages of /3-operators can be calculated using Wick-Bloch- 
de Dominicis theorem, e.g. 



^Kl, — K3 ^K2,—K4 ^Kl, 



s, 



K.2,—K,3 



1 + 1 + 1 + e^^-i 1 + 

/'kI /'k2'^K1, — K3'^K2, — K4 ~l~ f Kl f K2^ Kl , — K4^ K2 , — K.3 



(76) 



etc. After computation of these averages one finds that the coefficients 
near the averages contain the following products Xj^iJijj^ri-K) ^t+n,Kf^t-n^ 
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can be found with the help of (52). For simpUcity in what follows their values 
will be used in the case when J^^ — —J^^ — D. Then 



\+ ,,+ _ -inn _ \ - 



N 
1 



\+ ,,+ — p^Kn _ \ - - 

'^j+n,Kh*'j,—K T\T^ ^j+n,Kh*'j, — Ki 



\- ..+ 
jKr^j+rij — K 



N 
1 



-i(Kn—ipK.) 



N 



\- 11 + 



N 
1 



,i{Kn+tl}K.) 



(77) 



Gathering (74)-(77) together one derives the desired expression for transverse 
time-dependent correlation function for the model (30) 



4 < s'Ms 



j+n 



>- 



1 cosh {^—lEj, + inn + 

TV? cosh 



V 2 



1 sinh (^—tEi^t + tKn + iipi^ + 
TV? 



+ 



2 ; 



1 ^sinh(zV« + ^) 



H 2 



cosh (^) 



N 



X 



cosh (^) 



X 



1 sinh (y—iE^t + iKTi 
TV? 



cosh (^) 



(7J 



Although E^ and cos-^^, sin^^^ in (78) are determined by formulae (43), 
(63) for the case J^^ = — J^^ = D the obtained result covers the case (26) as 
well. Remembering formulae (27) and (31) one should simply use J^^cos^cc 
^j^y^y^ sin 2a + sin^ a instead of .F , J^^ sin^ a - i^^±^ sin 2a + 
JJ'?^ cos^ a instead of and i^l^^ instead of D with tan2a= ( J^?^ + J?^^) / 
I (^jxx _ jyyy jf Qj-^g puts D = in (78) it transforms into the well-known 
result obtained by Th.Niemeijer [21]. The depicted in figs. 17-20 depen- 
dence of the transverse dynamical autocorrelation function (78) on time 
{/3 = 10/ J^^; 1 : D = 0, n = 0, V : D = J^^, = 0; 2: 
D = 0, Vl = J^^, 2' : D = J^^, f2 = J^^) shows substantial changes 
caused by Dzyaloshinskii-Moriya interaction. 
The dynamical susceptibility 



Jo I ' ' 



(79) 
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is of great interest from the point of view of observable properties of the 
system. The obtained result (78) permits one to calculate the transverse dy- 
namical susceptibility. Really, taking into account the translation invariance 
one gets 

< [Sj{t), > = < ^|(^)4+n > - < Sji-tyj-n > = 



1 1 

— ^cos {ku - E^t) j^Y^sin [ku - Ei^t) tanh 



"2 



— ^cos {kti — E^t — i/j^) tanh 



+ 



-^^cos {ku - E^t + V'k) tanh '^^^^^ ^'^'^ ~ ^'^^ ~ ^'^^ 



1 



^sin {kti- E^t + il)^] 

K 



+ 



(80) 



you) 

Using for summation over sites in (79) the lattice sum Z^^=i e*'^" = 5^^, 
evaluating the integrals over t of the form 



cu + F^ + te' ^^^^ 

bearing in mind the definition of functions cosi/jk, sini/j^, and performing 
thermodynamical limit one ends up with 

. x_ 1 r , r 1 + 'p + ^ >-k) , 1 - cos (r, + r,_^) 
87r J-TT [Ep-K — Ep — uj — le —E^_p — Ep — uj — le 
1 + cos (V>p + -0p+^) _ 1 - cos (V>p + -0p+^) 
-Ep+^ + Ep-u;-ie E_p_^ + Ep - u - t 

Using the relation 



1 F 

tanh^. (82) 



= 'P-F^ + '""^^^P - 

Fp — ui — i£ Fp — u) 



(83) 



for real and imaginary parts of transverse susceptibility one gets final ex- 
pressions 

Rex..M = r dp [l+;^o^(^p + ^p--) + l- cos(^, + V.) _ 
87r J-TT [ Ep_^ — Ep — UJ —E^_p — Ep — ui 



1 + COS (jjp + ipp+^) _ 1 - cos (jjp + 1pp+^) 
-Ep+n + Ep-u E_p_^ + Ep-u 



tanh^, (84) 
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1 

Imxzz{>i, ^) = g y_ dp{[l + COS {jpp + V'p-k)] S{Ep-n -Ep-u) + 

+ [1 - cos (jjjp + V'p-k)] S(-E^_p -Ep-u)- 
- [1 + cos {ipp + V'p+k)] (5(--Ep+« + - u;) - 

- [1 - cos (^p + i/'p+«)] S{E_p_^ + Ep- w)} tanh 



^5) 



These are the main results of the present paper. 

It is useful to look at the particular case k — 0. In this case one has 



1 

Xzz{0,uj) = — J^^dpsin^ipp 



—Ep — E^p — u — IS 



Eo + E_ 



u — le 



tanh 



PEp 



(86) 



so that, for instance, 



1 

ImXzzi^,^) — ~~; / 0^9 sin^ '0p5(2£^p — a;) tanh 

4 J — TT 



PEp 



(87) 



In the case of isotropic XY model with Dzyaloshinskii-Moriya interaction 
simjjp = and /mx^^(0,a;) = as one should expect because in this case 
J2jLi Sj: H = 0. In the case of de Gennes model with Dzyaloshinskii- 
Moriya interaction when £p = \J^'^ + flJ cos p + one can integrate in 
(87) over p using the relation 



S{2£p - c) = ^ 



^(P-Po) 



Po 



2 








dp 





(88) 



where by po — po{u!) the solutions of the equation 2Spf^ —cu — are denoted. 
This equation can be written in the form 



cos Po = 

and when cu satisfies inequalities 



- 1< 



2 - J2 - 402 



< 1 



or for J > 



4nj 

\J -2n\<u; < J + 2fi, 



(89) 

(90) 

(91) 
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equation (89) has two sohitions in the interval of integration po ^ and 
— Po • Besides d£p/dp = —QJ sin p/2Sp, sinipp = J sin p/2Sp, so that in the 
case of de Gennes model with Dzyaloshinskii-Moriya interaction one gets the 
following final result 

(tanh ^ + tanh ^) , if| J - 2fi| < a; < J + 21], ^^^^ 
0, otherwise. 

The presented in figs. 21,22 results of the numerical calculations of frequency 
dependence of ImXzzi^,'-^) (86) for de Gennes model with Dzyaloshinskii- 
Moriya interaction (/3 = 10/ J^^; 1 : D = 0, 2 : D = 0.5J^^, 3: 
D — J^^) show that the presence of this interaction dramatically changes 
the frequency dependence. This fact seems to be of great importance in con- 
nection with the possible experimental prove of the presence in the system 
of Dzyaloshinskii-Moriya interaction on the base of experimental measure- 
ments of ImXzz{^,^)- 



6 Conclusions 

Let's sum up the results of present study of statistical mechanics of ID 
s = \ XY anisotropic ring in transverse field with Dzyaloshinskii-Moriya 
interaction. This interaction keeps the model in the class of ID s = g XY 
models because after fermionization of the Hamiltonian one is faced with the 
quadratic in Fermi operators forms. However, after their diagonalization one 
finds that the spectrum E,^ no longer is even function of k. This leads only 
to some technical complications in computations. The obtained thermody- 
namical functions and static spin correlation functions essentially depend 
on the value of Dzyaloshinskii-Moriya interaction. For instance, these in- 
teraction decreases the transverse magnetization at certain transverse field 
in de Gennes model and in isotropic XY model (figs. 5,6). They lead to 
appearance of non-zero spin correlators < s|sj_,_„ > and < sjsj_^„ > and 
to nonmonotonic dependence of pair spin correlation functions on n. The 
evaluation of transverse dynamical correlation function and the correspond- 
ing susceptibility shows that Dzyaloshinskii-Moriya interaction essentially 
influences on the dynamics of transverse spin correlations (figs. 17-20) and 
drastically changes the dynamical susceptibility (figs. 21,22). 

In addition it is necessary to note that if J^^ = J^^ = all obtained 
results transform into the corresponding results for anisotropic XY model. 
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Really, in this limit £k \J ei^^^ + 4(7++)^ sin^ k and E,^ = E^^- 

Due to this simplification the thermodynamical functions because of par- 
ity of integrands contain 2 Jq (k{. . .) instead of J^^ck{...). In contrac- 
tions (59)-(63) < (ffipf+n >^ Sn,0, < Vj^J+n >^ -<^n,0, < ^t^J+n >^ 

^ Jq (h COS {ku + ijjf^) tanh = G{n), < ^pjf'jj^^ >— > —G{—n) so that 
4 < sjs^_|_„ >, 4 < s^sj^n >, 4 < s|s|^,„ > (but not their squares) 
can be rewritten as x determinants and < „ >=< s^s?, „ >= 0. 
The transverse dynamical correlation function transforms into the obtained 
in ref. [21] expression. 

The performed investigations follow earlier works [18,22-29] considering 
the derivation of exact results in statistical mechanics of ID s — ^ systems 
with Dzyaloshinskii-Moriya interaction. 

At last it should be mentioned that for a quite a lot of magnetic and 
ferroelectric materials, showing nearly ID behavior above their ordering 
temeperatures, a variety of experimental data are now available [30-40] and 
thus theoretical investigations of statistical mechanics of ID spin models 
may be of great interest for clarifying whether the properties of such simple 
spin models are capable to caricature the measurements. 
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other participants of the Ukrainian- French Symposium "Condensed Matter: 
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